Equilibrium and off-equilibrium trap-size scaling in ID ultracold bosonic gases 
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We study some aspects of equilibrium and off equilibrium quantum dynamics of dilute bosonic 
gases in the presence of a trapping potential. We consider systems with a fixed number of particles 
and study their scaling behavior with increasing the trap size. 

We focus on one-dimensional (ID) bosonic systems, such as gases described by the Lieb-Liniger 
model and its Tonks-Girardeau limit of impenetrable bosons, and gases constrained in optical lattices 
as described by the Bose-Hubbard model. We study their quantum (zero-temperature) behavior 
at equilibrium and off equilibrium during the unitary time evolution arising from changes of the 
trapping potential, which may be instantaneous or described by a power-law time dependence, 
starting from the equilibrium ground state for an initial trap size. 

Renormalization-group scaling arguments, analytical and numerical calculations show that the 
trap-size dependence of the equilibrium and off-equilibrium dynamics can be cast in the form of 
a trap-size scaling in the low-density regime, characterized by universal power laws of the trap 
size, in dilute gases with repulsive contact interactions and lattice systems described by the Bose- 
Hubbard model. The scaling functions corresponding to several physically interesting observables 
are computed. 

Our results are of experimental relevance for systems of cold atomic gases trapped by tunable 
confining potentials. 

PACS numbers: 67.85.-d, 05.30.Jp, 67.85.Hj, 05.30.Rt 



I. INTRODUCTION 

The dynamics of strongly correlated quantum systems 
is a fundamental physical issue which has attracted much 
theoretical interest. The achievement of Bose-Einstein 
condensation in dilute atomic vapors [1] and the impres- 
sive progress in the experimental manipulation of cold 
atoms in optical lattices, see, e.g., Ref. Q and refer- 
ences therein, have provided a great opportunity to in- 
vestigate the interplay between quantum and statisti- 
cal behaviors in particle systems, and issues related to 
the unitary quantum evolution of closed many-body sys- 
tems, exploiting their low dissipation rate which allows 
to maintain phase coherence for a long time. Transi- 
tions between different quantum phases have been exper- 
imentally observed, such as those related to the forma- 
tion of a Bose-Einstein condensate in interacting Bose 
gases [l-Q and quantum Mott-insulator to superfluid 
transitions in atomic s yste ms constrained in optical lat- 
tices, see, e.g., Refs. [q-KLlj. Accurate experimental stud- 
ies of nonequilibrium properties of quantum many-body 
systems of ultracold atoms have also been reported, see, 
e.g., Refs. dMl- 

An important feature of these experiments is the pres- 
ence of a confining potential which traps the particles 
within a limited spatial region. The capability of vary- 
ing the confining potential, which may also depend on 
the spatial directions, allows to vary the effective spa- 
tial geometry of the particle systems, including quasi- ID 
geometries, see, e.g., Refs. |13l Il5l-ll9| . 

In the presence of an optical lattice with lattice spac- 
ing a, created by laser-induced standing waves which con- 
strains the particle to stay at the sites of a lattice, the the- 



oretical framework [2(| is provided by the Bose-Hubbard 
(BH) model [2l| with a confining potential V(r) coupled 
to the particle density, defined by the Hamiltonian 1 



(1) 



(ij) 



u 



where (ij) is the set of nearest-neighbor sites, hi are 
bosonic operators, rii = b\bi is the particle density op- 
erator, and N = (X)i n ») i s the particle number. We 
consider a power-law spatial dependence for the trapping 
potential, 



V(r) = -v p r p , 
P 



(2) 



where r = \x\is the distance from the center of the trap, v 
is a positive constant and p an even integer number. Ex- 
periments are usually set up with a harmonic potential, 
i.e., p = 2. Examples of experimental traps described by 
quartic potentials are reported in Ref. [22j. The hard- 
core (HC) limit U — > oo of the BH model implies that 
the particle number m per site is restricted to the val- 
ues rii = 0, 1. In one dimension the HC limit can be 



1 The BH Hamiltonian for N particles is usually written with the 
kinetic term "Hkin = — (J/ 2) Ulij) (bjbi + b\bj). The difference 
from Eq. JTJ is a ^-dependent constant. 
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exactly mapped into a lattice model of spinless fermions, 
see, e.g., Ref. [23j . 
We consider 



I = 



ji/p 



(3) 



as the size of the trap within the confined BH model [24 . 
[25j . In the case of harmonic traps, I ~ where uj is 
the trap frequency. The definition ([3]) of trap size nat- 
urally arises when we consider the thermodynamic limit, 
which is generally defined as N, I — > 00 keeping N/l d 
fixed and it is equivalent to introducing a chemical 

potential fi, adding the term /1 . n, to the Hamiltonian 
©• 

In the absence of a lattice structure, the basic model 
to describe the many-body features of a boson gas con- 
fined to an effective ID geometry is the Lieb-Liniger (LL) 
model with an effective two-particle repulsive contact in- 
teraction [23, 



u 



LL 



N 

E 

i=l 



2m 



V{ Xi 



i) (4) 



where N is the number of particles and V(x) is the 
confining potential. The limit of infinitely strong repul- 
sive interactions corresponds to a ID gas of impenetra- 
ble bosons HHH, the Tonks- Girardeau (TG) gas. ID 
Bose gases with repulsive two-particle short-ranged in- 
teractions become more and more nonideal with decreas- 
ing the particle density, acquiring fermion-like proper- 
ties, so that the ID gas of impenetrable bosons is ex- 
pected to provide an effective description of the low- 
density regime [3(J. ID Bose gases of cold atoms have 
been realized experimentally [lj, [H, GjHUi ■ 

In this paper we address issues related to the equi- 
librium and off equilibrium quantum dynamics of dilute 
bosonic gases in the presence of a confining potential. We 
consider ID trapped gases constituted by N bosonic par- 
ticles, at the equilibrium and off equilibrium during the 
unitary time evolution arising from changes of the trap- 
ping potential. In the latter case we consider a system 
of N particles in a trap of size Iq, which is prepared in 
its ground state at t = 0, the trapping potential is then 
varied as 



V(r,t) = -u p K(t/t q )r p , 
P 



(5) 



with K(0) = 1 and the parameter t q providing the time 
rate of the variation of the trapping. Since the particle 
number N = ^ x n x commutes with the time-dependent 
Hamiltonian even when the trapping potential depends 
on the time, the particle number N remains unchanged 
during the dynamical process. Istantaneous changes from 
the initial trap size Iq = 1 ju to a final trap size If, or 
a complete drop of the trap (corresponding to If — > 00), 
give rise to interesting cases of off-equilibrium evolutions. 



Moreover, we also consider a power-law time dependence 
of the confining potential such as 



K(t/t q ) = 



T = 1 + t/t q 



(6) 



Adiabatic time evolutions apply when the change of the 
external potential is very slow, thus it requires a large 
time rate \t q \. 

Several theoretical studies have already been dedicated 
to issues related to the quantum behavior of trapped 
bosonic gases, in the continuum and on the lattice (i.e., 
the BH model), in equilibrium conditions and off equi- 
librium due to time-variations of the trapping potential, 

see, e.g., Refs. gmMUJil- 

In this paper we further investigate these issues within 
the framework of the trap-size scaling (TSS) theory (251 
[59j . We consider systems with a fixed number of parti- 
cles N and study their scaling behavior with increasing 
the trap size I. We study the asymptotic trap-size depen- 
dence in the low-density regime which is characterized by 
universal power laws, as we shall see. 

Using renormalization-group (RG) scaling arguments, 
we derive the power-law behaviors which describe the 
asymptotic TSS in the low-density regime at equilibrium. 
They arc determined by the continuum nonrelativistic 
bosonic 4> 4 theory, which is the same continuum theory 
describing the quantum critical behavior at Mott transi- 
tions driven by the chemical potential [2lj . This implies 
that in one and two dimensions the power-law TSS is 
characterized by the dynamic exponent z = 2, the RG 
dimensions yb = d/2 and y n = d of the bosonic and 
particle density operators respectively [2ll m > an d the 
trap exponent 6 = p/(p + 2) [25]. This general TSS sce- 
nario is supported by analytical and numerical calcula- 
tions within ID systems. We show that the TSS is univer- 
sal for diluted gases with repulsive contact interactions 
such as the LL model (jlj and lattice systems described 
by the BH model (TTJ. 

We consider the off-equilibrium behavior arising from 
time variations of the confining potentials, which may 
be instantaneous or described by a power-law time de- 
pendence, such as Eq. ([6]), starting from the equilibrium 
ground state for an initial trap size Iq- We put forward 
scaling Ansatz for the asymptotic TSS with respect to 
the initial trap size lo in the large-Zo limit. Then we 
study the trap-size dependence of the off-equilibrium dy- 
namics of ID bosonic gases of N particles, by assuming 
an adabiatic evolution in the case of slow changes of the 
Hamiltonian parameters, and by analyzing the solutions 
of the Schrodinger equation of N impenetrable bosons in 
the presence of time-dependent harmonic potentials. 

Our results are of experimental relevance for systems 
of cold atomic gases trapped by tunable confining poten- 
tials. Indeed, the long characteristic time scales of these 
systems may allow a scaling study of the trap-size depen- 
dence of the zero-temperature properties of AT-particle 
boson gases in the low-density regime, at equilibrium and 
off equilibrium during the time evolution of the confining 
potential. 
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The paper is organized as follows. In Sec. [IT] we derive 
the asymptotic TSS of bosonic systems of N particles at 
equilibrium in the large trap-size limit, using general scal- 
ing arguments. In Sec. lIIII we focus on the ID BH model, 
we present analytical and numerical calculations in the 
HC limit and at finite U, which support the RG scaling 
predictions, and provide the universal scaling functions 
at equilibrium. In Sec. IIVI we consider ID bosonic gases 
at low density, and, in particular, in the TG limit which 
describes the low-density behavior of the LL model: we 
show that the trap-size dependence of TV-particle systems 
is identical to the asymptotic TSS of the ID BH model 
with the same number of particles. In Sec. [V] we turn 
to off-equilibrium dynamics, considering time-dependent 
confining potentials or instantaneous changes of its pa- 
rameters; using scaling arguments we extend the equilib- 
rium TSS Ansatz to the case of the off-equilibrium evolu- 
tions, distinguishing the cases of instantaneous variations 
and power-law time dependences of the confining poten- 
tial. In Sec. lVII we discuss the case of slow time variations 
of the trapping potential, assuming an adiabatic or quasi- 
adiabatic approximation. Sec. IVIII is dedicated to the 
study of the quantum unitary evolution of ID impenetra- 
ble bosons in time-dependent harmonic traps, where the 
off-equilibrium trap-size dependence can be analytically 
determined, in particular, for a linear time-dependence of 
the confining potential and for instantaneous quenches. 
Finally, in Sec. IVIlTl we summarize our main results and 
draw our conclusions. In App. [X] we discuss the case 
of particle systems confined by a spatial dependence of 
the hopping parameter. In App. |B]we study the asymp- 
totic behavior of the TSS functions for a large number 
of particles, showing that they have a nontrivial large- N 
power-law scaling. App. [C] presents a detailed analysis 
of the off-equilibrium dynamics of a quantum oscillator 
with a time-dependent frequency. 



II. TRAP-SIZE SCALING OF N PARTICLES AT 
EQUILIBRIUM 

Before discussing the TSS of N particles in the large-? 
limit, let us note that the large- 1 limit keeping N fixed 
differs from that performed at fixed chemical potential p, 
i.e., considering the BH Hamiltonian 

U lt = UBa + (p-J)^2n i . (7) 

i 

Indeed, the large trap-size limit, keeping p fixed, implies 
an increase of the particle number so that 

N/l d = p(fi/J) (8) 

asymptotically, where d is the spatial dimension and p(p>) 
is a finite function of p. This thermodynamic limit is 
usually considered when quantum transitions are studied 
in confined particle systems, see, e.g., Ref. f55j . 



We are interested in the low-density regime, which is 
related to the limit p — > p c where p{p c ) — 0, which corre- 
sponds to a low-density to empty-state transition, which 
may be considered as a n = Mott transition. In the ho- 
mogeneous BH model without trap, the low-energy prop- 
erties at Mott transitions driven by the chemical poten- 
tial p are described by a nonrelativistic U(l)-symmetric 
bosonic $ 4 field theory 21], whose partition function is 
given by 

Z = J [D(f>] exp ^- J ' dt d d x cj , 

£ = 0*9 i 0+-^|V0| 2 +r|0| 2 + U |0| 4 , (9) 
2m 

where r ~ p — p c . The upper critical dimension of this 
bosonic theory is d = 2, thus its critical behavior is of 
mean- field type for d > 2. For d = 2 the field theory 
is essentially free (apart from logarithmic corrections), 
thus the dynamic critical exponent is z — 2 and the RG 
dimension dimension of the coupling p is = 2. In 
d = 1 the theory turns out to be equivalent to a free field 
theory of nonrelativistic spinless fermions, thus z = 2 
and Up = 2 as well, see, e.g., Ref. [23[ 

The quantum critical behaviors in the presence of the 
trapping potential can be described in the theoretical 
framework of the TSS theory (2||, which introduces a 
trap critical exponent 9 which describes how the length 
scale at the quantum critical point diverges with increas- 
ing the trap size I, i.e., £ ~ l e [Hj]. The trap exponent at 
the Mott transitions of ID and 2D BH models is 

e = P /(p + 2). (io) 

At the low-density to empty-state transition the TSS of 
the free-energy density in the presence of a confining po- 
tential © is given by [H [55| 

F(p, T, /, x) = r 9{ - d+z) F{pl e/u , Tl 9 \ xr e ), (11) 

where x is the distance from the middle of the trap, T 
is the temperature, p = p — p c , and v = l/y M . The 
zero-temperature TSS of a generic observable, whose low- 
density critical behavior is described by the RG dimen- 
sion y a in the homogeneous system, is given by 

(0){p,l,x) - l-Vo e O{jll 6 ' v ,xl- e ). (12) 

For example, any low-energy scale at T — is expected 
to behave as E = l~ ze £{pl 9 / v ); the particle density as 
(n x ) = l- de V(pl e ^,xr e ) (using the fact that the RG 
dimension of the density operator is y n — d); the one- 
particle density matrix as 

Pl (x,y) - (bib y ) - l- de M{pl e '\xl-\yl- e ), (13) 

using the fact that yt — d/2; etc.... 

The limit p — ¥ corresponds to the low-density regime 
Na d /l d -> where a is the lattice spacing. A TSS Ansatz 



4 



for the large-/ trap-size dependence at fixed particle num- 
ber N can be derived by replacing the dependence on 
\±l B l v with that on N. See also the next section for an 
explicit derivation in ID systems. Therefore, for a generic 
observable we expect 

(0)(N, Z, x) ~ r y ° e O N (xr e ). (14) 

In particular, the gap, i.e., the energy difference of the 
lowest states, behaves as 

A N « A N l-' e , (15) 

where An is a ./V-dependent amplitude. Since the RG 
dimension of the particle density operator rii is given by 
y n = d + z — y M = d, we expect that 

p(x) = (n x ) « r de V N (X), (16) 

and 

G n (x,y) = (n x n y ) - (n x )(n y ) a l~ 2d6 Q N (X, Y), (17) 

where X = x/l e and Y = y/l e , and the scaling functions 
T>m and Gn depend on N. The RG dimension of the 
boson operator is y\, — d/2, thus the low-density TSS of 
the one-particle density matrix is 

Pl (x,y) = (btb v ) « r de M N (X,Y). (18) 

Finally, we consider the momentum distribution, defined 
as 

normalized so that J ^rik = 1. tt.^ is usually accessi- 
ble experimentally from the intereference patterns of ab- 
sorption images taken after the drop of the trap and the 
expansion of the atomic gas. Eq. (IT8l) implies 

n k ^l e{2 - d) Af N (K), K = l 9 k. (20) 

In the next sections we report analytic and numerical 
calculations for the BH model and a gas of impenetrable 
bosons, showing that they share the same asymptotic 
TSS behavior, with the same scaling functions. 

III. TRAPPED PARTICLES AT EQUILIBRIUM 
ON A ID LATTICE 

In this section we address issues related to the equi- 
librium properties of a ID lattice system of TV bosonic 
particles confined by a trapping potential, described by 
the ID BH model |T]). These results will also be rele- 
vant for the off-equilibrium evolution of the system under 
variations of the confining potential, in particular when 
the dynamics is so slow to admit the adiabatic approxi- 
mation. In the following we set the hopping parameter 
J = 1, thus the trap size Q simply becomes I — 1/v. 



To complete the definition of the BH model in the pres- 
ence of the confining potential, we consider traps whose 
center coincides with a site of the lattice, so that the lat- 
tice model has a reflection symmetry with respect to the 
center of the trap. However, any other particular choice, 
i.e., centering the trap anywhere between two sites, does 
not change the asymptotic low-density TSS behavior, 
leading to an effective asymptotic reflection symmetry. 

A. The hard-core limit 

We first consider the hard-core (HC) limit U — > oo of 
the ID BH model, which allows us to study the effects 
of the confining potential by exact and very accurate nu- 
merical results. The HC limit implies that the particle 
number per site is restricted to the values rii = 0,1. 
In this limit the ID BH model ([7]) can be mapped into 
the XX chain model with lattice spacing a and a space- 
dependent transverse external field, 

H X x = - ^ {SfSf +1 + SfSf +1 ) 

i 

- £> + V(*i)]Sf, (21) 

i 

where Sf — of/2 and a a are the Pauli matrices, which 
are related to the boson operators bi by of = b] + 6,, 
of = i(b\ — bi), of = 1 — 2b\bi. Then, by a Jordan- Wigner 
transformation, one can further map it into a model of 
spinlcss fermions, see, e.g., Ref. [23|, with a quadratic 
Hamiltonian 

H q = J2 c ihiiCj + (M - 1) E C K (22) 

ij i 

where a is a spinless fermion operator, and 

hij = Sij - 7^;j-i - 7^ij+i + V(x l )5 ij . (23) 

Issues related to quantum transitions in particle sys- 
tems are best discussed in the presence of the chemical 
potential /i. In the absence of the trap, the ID HC-BH 
model with a chemical potential /i has three phases: two 
Mott insulator phases, for /i < — 1 with (m) — 1 and 
for fi > 1 with (rii) — 0, separated by a gapless super- 
fluid phase for \fi\ < 1. Therefore, there are two Mott 
insulator to superfluid transitions at fj, = — 1 and fi = 1. 
These transitions are characterized by the dynamic ex- 
ponent z = 2 and the RG dimension of the chemical 
potential j/ M = 2 [2l|. The gapless superfluid phase is 
instead described by a free massless bosonic field theory 
with dynamic exponent z = 1, see, e.g., Ref. [231 ] . 

In ID particle systems, the thermodynamic limit at 
fixed /i corresponds to N, I — > oo keeping the ratio N/l 
fixed. Indeed, we have 

N = (J2bth)=p(fx)l + 0(l) (24) 
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The function p{p) can be computed in the HC limit. The 
particle density in the large-^ limit turns out to approach 
its local density approximation (LDA), with corrections 
that are suppressed by powers of the trap size and present 
a nontrivial TSS behaviour 55]. Within the LDA, the 
particle density at the spatial coordinate x equals the 
particle density of the homogeneous system at the effec- 
tive chemical potential 

fl eS (x) =fM+- (y)*\ (25) 

The LDA of the particle density reads (n x )n a = 
Pida(z/Z), where 

{0 for fies(x) > 1, 

(1/tt) arccos^ c ff(x) for -1 < p, eS (x) < 1, 
1 for fiee(x) < -1. 

(26) 

Asymptotically, the total particle number is obtained by 
integrating the LDA of the particle density pida, obtain- 
ing 

pifi) = 2 / p Ua {y) dy. (27) 
Jo 

In the low-density regime, p, = p, — 1 — > 0, 

^)=c|/2|( 2 +p)/^)[1 + 0(m)] ) (28) 

where c is a p-dependent constant; c = 1 for p = 2. 

Eq. (j2"4")l provides the correspondence between the ratio 
iV/Z and p. In particular, when < N/l < p(—l) the 
system is effectively in the superfiuid phase, while for 
N/l > p{— 1) the n = 1 Mott phase appears around the 
center of the trap. Eq. fl27J gives /5(-l) = 2.54648 for 
p = 2 and p(-l) = 2.56561 for p = 4. Eqs. (JUJ), CT . 
and ([2§| imply that the particle density at the origin 
scales as a nontrivial power law of the ratio N/l in the 
low- density regime, 

p(0) ~ (AT/Z)*. (29) 

At the low-density Mott transition, around p, = 1, 
the TSS limit can be analytically derived within the 
quadratic spinless fermion representation (55| . This is 
obtained by rescaling the spatial distance from the ori- 
gin as 

x = l e X, (30) 
and the difference p, = 1 — \x as 

£ = r 2 V- ( 31 ) 

Any low-energy scale turns out to behave as E rj 
l~ 2e £A{Hr), the particle density behaves as (n x ) w 
l~ e T>(p r , X), etc., in agreement with the scaling Ansatz 
reported in the previous section, cf. Eq. (jT2"j) . Analytic 



and numerical calculations of the above scaling functions 
are reported in Ref. [55j. 2 

B. TSS of a system of N particles 

We now derive the TSS as a function of the particle 
number N in the low density regime Na/l <C 1. This 
is worth being discussed in some detail, because the cor- 
responding TSS functions are not trivially derived from 
the TSS in the presence of a chemical potential, which 
were already computed in Ref. [55[. 

1. TSS limit 
In the fermion representation the Hamiltonian 

»><••'•.• ( 32 ) 

ij 

can be diagonalized by introducing new canonical 
fermionic variables r\k = "^2i4>kiCi, where <f> satisfies the 
equation 

hij4>kj = u k 4> k i, (33) 

so that 

H c = y^^kvtvk- (34) 

k 

The ground state of a system of N particles is then given 
by the ry-fermions filling the N lowest one-particle lev- 
els. For the lowest states, assuming smoothness, we may 
consider the continuum limit of Eq. (|3"3"|) . by replacing 
4>kx — > <f>k(x) and rewriting the discrete differences as 

, . , , . . dd>(x) 1 „d 2 6{x) , „ 

<t>(x + a)-<f>(x) =aJ ^ i + 2 a + ( 35 ) 

where a is the lattice spacing. Then, rewriting the re- 
sulting equation in terms of the rescaled quantities 

X = a- 2e/p r e x, (36) 
e k = a- 2e l 2e uj k7 (37) 
Vk {X)^a s ln e l 2 ^ k {a- 2S ln e XX (38) 

with 8 given by Eq. fj 10[) , and neglecting terms which are 
suppressed by powers of the trap size, one arrives at a 
Schrodingcr-likc equation 

(-l^ + ^xA<p k (X) = e Wk (X). (39) 



2 Note that here the definition of the trap size, cf. Eq. (2), differs 
from that of Refs. by a factor p" 1 ^. 
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This equation describes the TSS limit at fixed N, i.e., 
I — > oo, x — > oo, keeping the scaling variable X fixed. The 
next-to-leading terms in the large-Z limit, arising from 
the higher order terms in the expansion (|35l) . give rise 
to 0{l~ 20 ) scaling corrections [55(. Moreover, one can 
easily check that a shift of the center of the trap, by 
5 < a, generally induces 0(l~ e ) subleading corrections. 
Solving Eq. (|39|) for p = 2, we obtain 



e k = k + l/2, fc>0, 
fk{X) - 



(40) 



7r i/42Va(jfe!)V5 



exp(-X 2 /2), 



where X = xf(al) 1 ^ 2 and are Hermite's polynomi- 
als. For p = 4, Eq. can be solved numerically 
by Numerov's method, see, e.g., Ref. [gfj; the result- 
ing energy levels are e = 0.420805, e x = 1.50790, 
e 2 = 2.95880, e 3 = 4.62122, e 4 = 6.45350, e 5 = 8.42843, 
etc.... The Bohr-Sommerfield quantization formula, see, 
e.g., Ref. [6lj], gives the asymptotic large- k behavior 



h(k + 1/2) 4 / 3 , 

7r 2 / 3 r(7/4) 4 / 3 

2 4 /3r(5/4) 4 /3 



0.867145. 



(41) 



This formula provides a good approximation for rela- 
tively low levels already: it is accurate to 0.1% already 
for e$. For p —> 00, Eq. (|39| becomes equivalent to the 
Schrodinger equation of a free particle in a box of size 
L = 21 with boundary conditions tp(— 1) = <£>(1) = 0, 
leading to 



fc > 0, 



(42) 



^(X) = sin -{k + l)(X + l) 



where X = x/l. 

In App.[3]we show that the same TSS limit is obtained 
when the trap is induced by a spatial dependence of the 
hopping parameter. 



2. TSS of observables 

In the following we set the lattice spacing a = 1 to 
simplify the expressions. The dependence on a and J 
can be easily recovered by a dimensional analysis. 

The TSS limit leading to Eq. (1591 allows us to compute 
the low-density trap-size dependence of the observables. 
For example, the gap, i.e., the difference of the energy of 
the lowest states, behaves as 



A N = A N r™ [i + o(r a % 

with A2 = 1 for p = 2, and 

A N = apN 20 - 1 [1 + 0(1/N)} 



(43) 



2 
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FIG. 1: (Color online) The scaling function Vn(X), cf. 
Eq. (HI, for p = 2. Since X>jv(X) = £>iv(-X) due to the 
reflection symmetry with respect to the center of the trap, we 
show only the curves for X > 0. 
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FIG. 2: (Color online) The scaling function g N (0,X), cf. 
Eq. (JUJ), for p = 2. 



for a generic power law p: 04 1.15619 and = ir 2 /2. 
Note that the gap at fixed N differs from the difference 
of the energy of the lowest states at fixed chemical poten- 
tial n, which behaves as A M = l~ 2 £a(/v)i because the 
latter involves states of subsequent particle number sec- 
tors, giving rise to zeroes in the scaling function £a(m»0 
for [i r < 0, see Ref. (55j . 

In the low-density regime the particle density behaves 

as 

p{x) = (n x ) = r e V N {X) [1 + 0(r 26 )} , (45) 

N-l 

V N {X) = J2 ti{X). (46) 

Fig. \H shows results for the spatial dependence of 
Vn(X) for p = 2 and several values of N. Note the 
peculiar structure of T>n(X) characterized by N local 
maxima, which get suppressed at large N by powers of 
l/N. Due to the parity of the Hermite polynomials, 



(44) V 2j _ 1 (0)=V 2j (0). 
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FIG. 3: (Color online) Results for the one-particle density 
matrix: we plot I pi(0,x) vs X = x/l e for p — 2, thus = 
1/2, for several values of N. The data points are numerical 
results for the HC-BH model at fixed N and trap size with 
10 < I < 10 3 . The continuous lines are the curves for systems 
of N impenetrable bosons. The data of the HC-BH model 
clearly approach these curves in the large trap-size limit. 



Straightforward calculations show that the density- 
density correlator behaves as 

G n (x, y) = (n x n y ) c « r 2e G N (X, Y), (47) 

where X = x/l e , Y = y/l e , and 



QN{X,Y) = -\E*£ip k {X)<p k {Y)\ . (48) 

Fig.[2]shows plots of Qjsr(0, X) for the harmonic potential. 

The one-particle density matrix, cf. Eq. (fi"8|) . cannot 
be easily derived from the solutions of Eq. , because 
the fermion-boson map exploited in the HC limit is not 
trivial, and, in particular, it is non local. However, as 
we shall show in Sec. IIV1 the asymptotic TSS of the 
BH model in the low-density regime coincides with the 
trap-size dependence of a ID gas of impenetrable bosons, 
whose one-particle density matrix can be computed us- 
ing the known ground-state wave function. Some results 
for the harmonic potential are shown in Fig. [3J 

The TSS functions of the observables considered above 
show nontrivial power-law scalings with respect to the 
particle number N at large N. Their large- N behaviors 
are reported in App. |Bj 



C. Numerical results 

Beside deriving the asymptotic behaviors in the low- 
density region, we present numerical calculations at fixed 
particle number N and trap size I. We exploit the 
quadratic spinless fermion representation (l22t of the ID 
HC-BH model, which allows us to perform computations 
for very large systems, since they only require the diag- 
onalization of a L x L matrix where L is the number 



FIG. 4: (Color online) IA N versus N/l for p = 2. The dashed 
horizontal line indicates the constant value computed in the 
low-density regime, i.e., N/l -C 1. The vertical dotted line 
shows the asymptotic value of the ratio N/l corresponding to 
the n = 1 Mott transition. 
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FIG. 5: (Color online) Plot of N~ irA l A/3 l^ N versus (iV//) 4/3 
for p — 4. The dashed horizontal line indicates the constant 
value computed in the low-density regime, i.e., N/l <S 1. The 
vertical dotted line shows the asymptotic value of (N/l) i ^ 3 
corresponding to the n = 1 Mott transition. 



of lattice sites. We obtain numerical results for chains of 
size L, with a trap of size / centered at the middle site (we 
consider odd L); we choose L large enough to have neg- 
ligible finite-L effects. This can been accurately checked 
by comparing results at fixed I and increasing values of L. 
Thus, the results at fixed N and I that we shall present, 
respectively up to N w 10 2 and I — O(10 3 ), are the in- 
finite chain size limit (keeping N and I fixed) with great 
accuracy. For more details see Ref. (55[ , where analogous 
calculations at fixed chemical potential were presented. 

Bosonic particle systems confined to ID lattices have 
already been the subject of several numerical investiga- 
tions H3, IH-Sl |5l[54j. We study the dependence 
of some physically interesting observables on the particle 
number N and the trap size I. 

Figs. Q] and [S] show results of the gap for p = 2 and 
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FIG. 6: (Color online) The particle density at the origin ver- 
sus (N/l) 1 / 2 , for p — 2 for several values of N and I. The data 
points show results obtained by solving the HC-BH model at 
fixed iV and I, while the continuous line shows the LDA. The 
vertical dotted line shows the asymptotic value of (N/l) 1 ^ 2 
corresponding to the n = 1 Mott transition. 
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FIG. 7: (Color online) The particle density at the origin 
versus (N/l) 2 ^ 3 , for p = 4. The data points show results 
obtained by solving the HC-BH model at fixed and I, while 
the continuous line shows the LDA. The vertical dotted line 
shows the asymptotic value of (N/l) 2 ^ 3 corresponding to the 
n — 1 Mott transition. 



p = 4 respectively. Note that the gap of a system of 
TV impenetrable bosons is identical to the gap of N free 
fermion particles in a trap. Guided by the low-density 
scaling behavior (|43[). we plot the quantity N 1 ~ 2e l 26 IS.^ 
versus (N/l) 26 (which is just IA N vs. N/l for p = 2). 

In the region of values of N/l corresponding to 1 > 
fi > — 1, the data show the asymptotic behavior 



Around N/l = p(—l), i.e., the value corresponding to 
the n — 1 Mott transition, the behavior for even and odd 
N begins differing significantly. In particular, the data 
for even N appear suppressed for N/l > p(—l). This is 
essentially related to the fact that the trap is centered at 
the middle site of the chain. When the region around the 
center of the trap shows the n = 1 Mott phase, we have 
two degenerate lowest states for even N, differing for a 
reflection with respect to the middle site; while for odd 
N the ground state is unique, and the gap is expected 
to behave as A at ~ NP^ 1 /l p for N/l sufficiently larger 
than p{— 1), as also shown by Fig. 2]for p — 2, where the 
corresponding asymptotic behavior I An ~ N/l can be 
already seen for N/l > 3. 

Figs. |B]and [7] show results of the particle density at the 
origin, for p — 2 and p — 4 respectively, for some values 
of N in the range 20 < N < 100, and I up to O(10 3 ). 
The data appear to follow a unique function of p = N/l, 
given by the LDA obtained using Eqs. (|2"6) and (|2"T1) . and 
corrections are hardly visible. This fact was already ob- 
served by other numerical works, see, e.g., Refs. PHI. 



The agreement is already good for relatively small values 
of N, i.e., N > 20. Note that the LDA reproduces the 
low-density behavior (|45|) for N ^> 1, and in particular 
the leading large- N term of Eq. (IBip . 

Figures |S] and [3] show results for the momentum distri- 
bution Ufc, cf. Eq. (fll?l) . for p — 2 and p = 4 respectively, 
and several values of N and I. We first note that the 
plots of nk vs k show the scaling behavior 



n k ^f(N/l,k) 



(51) 



Moreover, as shown by the plots of l~ e rik versus K = l e k, 
the data also support the TSS behavior (|20|) . which is 
expected to be approached with 0[(N/l) ] corrections. 
Actually, as shown by the bottom figures © and ©, the 
data for k > appear to scale as 



(52) 



n k - {N/l)-°F{K), K ee (N/l)~ e k, 



for N/l < 1, which agrees with both Eq. ^STJ) and ([20 
The zero component 



n 



(53) 



scales differently, indeed no = 0(1) in the large- N limit, 
analogously to a gas of impenetrable bosons. See, e.g., 
Ref. [Hj]. At large K, F(K) - w hi c h can be in- 

ferred from the results of Refs. [H^, Ej] for a gas of im- 
penetrable bosons. 



NA N = g(N/l), 



(49) 



in the limit / — > oo, AT — > oo keeping AT/i fixed. The low- 
density behavior (|4U)) is recovered for A^/Z <C 1, because 



cx 



20 



(50) 



D. Finite U and universality of the low-density 
behavior 

We now consider the BH model at finite values of the 
on-site repulsion coupling U. We perform calculations 
using the DMRG method. Specifically, we consider the 
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FIG. 8: (Color online) Results for the momentum distribution 
rife for p — 2, and several values of N and I. We plot nt vs 
k (above), l~ e nk versus K = l B k (middle), and (N/l) B rik vs 
(N/l)~ e k (below). We recall that = 1/2 for p = 2. 
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FIG. 9: (Color online) Results for the momentum distribu- 
tion rife for p — 4, and several values of N and I. We plot rife 
vs k (above), l~ e rik versus K = l B k (middle), and (N/l) 7ik 
vs k/(Nl) e (below). We recall that 6 = 2/3 for p = 4. 



BH model with U = 2 in the presence of a harmonic 
potential, up to trap sizes I — O(10 3 ), and for several 
values of N, up to N — 20. The trap is again centered 
in the middle site of a lattice of size L, which is taken 
sufficiently large to make finite-L effects negligible. We 
set the cutoff on the number of bosonic states per site 
n b = 5, which turns out to be sufficient to provide very 
accurate results; indeed the relative difference from the 
results using ub = 4 is at most O(10~ r ). We keep the 
maximum eigenvalue truncated from the density matrix 



below 10 -10 ; this requires retaining up to 200 states. 

The issue that we want to investigate is the univer- 
sality of the low-density TSS behavior with respect to 
variations of the on-site repulsion coupling U, i.e., how 
it depends on U. If there is universality, then the low- 
density asymptotic TSS at finite values of U must be the 
same as that found in the HC limit. Actually, a reseat- 
ing of the trap size which depends on U may be allowed, 
although the data show that this is not necessary when 
one uses the same definition of trap size with respect to 
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FIG. 10: (Color online) Some results for the gap of the BH 
model with (7 = 2 and TV = 5, 6, 10, 20. They show that the 
data approach the asymptotic value I An = 1 with increasing 
I, with 0(/ _1 ) corrections. 
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FIG. 11: (Color online) The difference l 1/2 {n ) - ©jv(O) vs 
{N/l) 1/2 , for the BH model with U = 2 and TV = 5, 6, 10, 20. 
These results provide a clear evidence that the asymptotic 
value is consistent with Djv(0), and it is approached with 
0(i~ 1//2 ) corrections. 



the kinetic term, as we have done in Eqs. ([JJ and ©. 

Fig. [10] shows results for the gap, i.e., the difference 
between the energy of the lowest states. The data show 
a behavior analogous to that found analytically in the 
HC limit, see Sec. IIIIB21 i.e., IA N = 1 + 0{l~ x ) inde- 
pendently of the particle number 7Y [the large-Z extrap- 
olation to get the leading behavior is checked within an 
accuracy of O(10 -6 )]. This shows that there is no need 
of a {/-dependent normalization of the trap size. 

Fig.lTTlshows data for the particle density at the origin. 
They are consistent with 



l 1 / 2 p(0) = V N (0) l + oii- 1 / 2 ) 



(54) 



where T> N (0) = \/2N/ir, using Eqs. flU]) and (gDJ). Note 
that the power law of the scaling corrections differs from 
that found in the HC limit, which was 0(l~ 26 ), thus 
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FIG. 12: (Color online) Some results for the spatial depen- 
cence of the particle density in the BH model with U — 2 
and for TV = 5 (above) and TV = 10 (below). They clearly 
approach the large-/ limit of the HC model, although their 
convergence appears significantly slower than that found in 
the HC limit. 



0(l~ r ) for p = 2. Therefore the approach to the asymp- 
totic behavior is significantly slower than that for the HC 
limit. This is also found for other observables. Results for 
the spatial dependence of the particle density are shown 
in Fig. [T^lfor TV = 5 and TV = 10. They clearly approach 
the scaling function T>n(X) with increasing the trap size. 
Finally, in Fig.[T3]we show results for the momentum dis- 
tribution at several values of TV, and compare them with 
the asymptotic behavior computed in the HC limit. 

In conclusion, these results confirm that the low- 
density TSS of TV particles described by the BH model is 
universal with respect to the on-site repulsion coupling 
U. However, scaling corrections at finite U appear gener- 
ally larger than those of the HC limit, 0(l~ 6 ) for generic 
values of U against 0(l~ 28 ) in the U — > oo HC limit. The 
only exception was the gap where we have not found ev- 
idence of 0(Z -1 / 2 ) corrections. 

Actually, 0(l~ e ) corrections are generally expected, 
because the RG dimension of the parameter U has RG 
dimension yu = — 1 with respect to the low-density scal- 
ing. This implies that it generally leads to 0(^ _1 ) scaling 
corrections, which become 0(1 ) in terms of the trap 



11 



10 



10 



= 20 
= 50 
= 100 
= 250 
-- 500 
= 1000 
= 2000 



HC 



I n 




FIG. 13: (Color online) Some results for the BH model with 
U = 2 for TV = 5 (above) and N = 10 (below). 



size. 3 These corrections vanish in the HC limit. Thus, 
within generic ID BH models, the HC limit represents a 
RG-improved model [gj] where the leading scaling cor- 
rections are absent. 



IV. THE ID BOSONIC GAS AT LOW DENSITY 

We now consider a system of ID boson particles inter- 
acting through a repulsive contact term in the presence 
of a confining potential such as ©, described by the LL 
model, cf. Eq. @. In the low-density regime the sys- 
tem can be effectively described by the limit of infinitely 
strong repulsive interaction |3G |. i.e., a ID gas of impen- 
etrable bosons (TG model). 

The wave function for a ID system of N impenetrable 
bosons in a confining potential is essentially defined by 



3 The RG dimension of U can be derived from the /3-function 
associated with the quartic term of the corresponding bosonic 
continuum theory which reads 0(u) = u — u 2 /2 exactly :23j. 
It has a nontrivial fixed point for u* = 2, thus y\j = /?'(«*) = —1. 



the one-particle Hamiltonian 

» 2 1 
H = + -rwJ > x v 
2m p 



(55) 



and the impenetrability condition, i.e., the fact that the 
wave function of the N particles vanishes if two spa- 
tial variables coincide. The low-density condition to 
realize a ID TG gas of impenetrable bosons is @, H(| 
Na 2 s /l 2 sc << 1 where a s is the ID scattering length, re- 
lated to the quartic coupling by g = —AH 2 /(m 2 a s ), and 



In 



«Va/( 



,1/2 



is the oscillator length. In the fol- 
lowing we set h — 1 and m = 1. 

The ground-state wave function of the TG model can 
be written in terms of the ground state wave function of 
N free fermion particles [HI, [29[ described by the Hamil- 
tonian (l55l. which is 



...,x N ) 



1 



fN\ 



(56) 



where 4>i(x) are the lowest N eigensolutions of the one- 
particle Schrddinger equation Hfy — E^. The wave 
function $ of N impenetrable bosons is obtained by sym- 
metrizing the fermion wave function i.e., 

${x\,...,x N ) = A(xi,...,x N )^{xx,...,x N )-, (57) 
A(xx,...,x N )= Yl siga(xi-Xj). (58) 

l<i<j<N 

The ground-state wave function allows us to derive the 
the one- and two-particle density matrices by 



Pi(x,y) = N I $(x,x 2 ,...,x N )*$(y,x 2 , ...,x N )dx 2 ...dx N 

(59) 



and 

p 2 {xi,x 2 ;y 1 ,y 2 ) 



(60) 



N 



$(xi,x 2 ,x 3 , ...,x N )*$(yi,y 2 ,x 3 , ...,x N )dx 3 ...dx N . 



In the case of the harmonic potential, we have 

E k =w(fc + l/2), fc>0, (61) 
<PkW u nl/42k/2{M)1/2 e 

thus leading to [35| 

■5>{x ll ...,x N ) = c N uo N2 l i B(x 1 ,...,x N )e-^^l 2 , 



B(xi,...,x N ) 



n \ x i~ x ii 

\<i<j<N 



(62) 



where cjv is the appropriate normalization constant 

n -1/2 



C N =7f 



-N/4: 



N-l 



fe=0 



(63) 
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so that J Y\ i=1 d,Xi\&\ 2 = 1. Some useful analytical de- 
velopments, to evaluate the one-particle density matrix, 
can be found in Refs. [H EH EI EI EI H3] ■ 

We now note that, after appropriate rescalings, the 
above results for the TG model concide with the low- 
density TSS of the BH model, see Sec. IIII B 1[ which was 
derived by taking the continuum TSS limit of its HC 
limit. This implies that the trap-size dependence of the 
TG model exactly gives the asymptotic TSS of BH model, 
after replacing 4 I = a; -1 . This can be verified by explicit 
calculations, see below. 

Straightforward calculations lead to the following ex- 
pressions for the particle density and its correlator: 

p(x) =px{x,x) =r e V N {X) (65) 

and 

G n (x,y) = (n x n y ) c = (66) 
= Pi {x, y, x, y) - pi (x, x)px (y, y) 

= r 29 g N (x,Y) 

where 9 = p/(p + 2) is the trap exponent already intro- 
duced in the low-density TSS of the BH model, X = x/l e 
and Y = y/l e , and the TSS functions T>n and Gn are ex- 
actly given by Eqs. (|46|) and (|48)) . 

The one-particle density matrix can be written as 

Pl (x,y) =r e M N {X,Y). (67) 

Again, this exactly provides the large-^ TSS of the BH 
model, as shown by results for the TG model and 
the HC-BH model in Fig. [3l In particular, its large- 
N limit is given by Eq. (|B6j) . which implies that the 
rescaled density matrix (N/l)~ e pi(x, y) has a nontriv- 
ial large- N limit B((,S) keeping £ = xN- 1+e l^ e and 
S = (y — x)N~ e l~ e fixed. This scaling behavior was al- 
ready noted in Ref. [42l |. 

We can also compute the energy difference A^r be- 
tween the two lowest states. The lowest excited state 
above the ground state is obtained by exciting only the 
fermion particle with the highest energy in the ground 
state. One can easily check that A at is exactly given 
by the asymptotic TSS behavior found for the HC-BH 
model, cf. Eq. (I43j) . without corrections. In particular, 
A^v = 1/1 for the harmonic potential. 



4 Restoring the dependences on J, a and m in the BH and TG 
models, the trap-size correspondence between the trap sizes of 
the BH and TG models is 



thus I = ui~ 1 setting J,a,m to one. Therefore, the trap size 
of the TG model in a harmonic potential is essentially given 
by I ~ %/{mu), thus I ~ 1% SC where i osc = ^/h/(mui) is the 
characteristic length scale of an oscillator of frequency uj. 



Summarizing, we have shown that the trap-size depen- 
dence in a ID trapped gas of N impenetrable bosons coin- 
cides with the asymptotic TSS of N particles described 
by the ID BH model, if appropriate definitions of the 
trap size are considered. As already argued within the 
BH model, the critical exponents associated with this 
TSS are related to the nonrelativistic $ 4 theory We 
expect that the low-density TSS is also universal with re- 
spect to the strength of the short-ranged repulsive inter- 
action. Therefore, it should exactly provide the asymp- 
totic low-density trap-size dependence of N boson parti- 
cles described by the LL model, when Na 2 /l 2 sc <C 1. 

Note that the power-law TSS does not have corrections 
in trapped systems of impenetrable bosons, while within 
the BH model it is only expected asymptotically in the 
large-Z limit, i.e., it is approached with O{1~ 20 ) correc- 
tions in the HC limit and 0(l~~ e ) corrections for finite U. 
In a sense, in the language of the RG theory [65j|, the 
TG model represents a fixed-point Hamiltonian, i.e., a 
model where scaling corrections are totally absent, with 
respect to the low-density behavior of the BH model and 
the LL gas. Using the same RG arguments reported at 
the end of Sec. IIII Dl we predict that scaling corrections 
are 0{l~ e ) in the LL model. 



V. TRAP-SIZE SCALING IN A 
TIME-DEPENDENT TRAP 

The off-equilibrium dynamics is a quite complicated 
issue, more subtle than issues related to the equilibrium 
behavior. This is not a prerogative of the quantum evo- 
lution only, but it is also found in classical systems, see, 
e.g., Refs. [iE H3l- 

In this section we discuss the trap-size dependence of 
the off-equilibrium time evolution of ID bosonic gases 
in time-dependent traps, in the limit of instantaneous 
variations and for a power-law time dependence, starting 
from the equilibrium ground state for a initial trap size 
^o- We derive scaling Ansatz for the asymptotic TSS with 
respect to the initial trap size Iq in the large-Zp limit. 



A. Instantaneous variation of the confining 
potential 

Let us first discuss the case of an instantaneous change 
of the confining potential. In particular, we assume that 
at t = the TV-particle system is at equilibrium, in the 
ground state with a confining potential of trap size lo. 
Then, the trap is instantaneously changed to a larger 
trap size, If > Iq, or dropped completely, corresponding 
to // = oo. We are interested in the asymptotic trap- 
size dependence of the quantum time evolution after the 
instantaneous quench for large initial trap size lo. 

In this case we expect that the trap-size dependence 
of the off-equilibrium dynamics after the quench is es- 
sentially determined by the trap-size dependence of the 
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inizial state at equilibrium, and by the ratio lf/lo of the 
final and initial trap sizes. Thus, the simplest scaling 
Ansatz for the large-/o behavior may be 



(0) N (x;t) a tf" e O N (xlo e ,tlo M0 ) lf/lo), 



(68) 



where z = 2 and 9 = p/(p + 2) is the equilibrium trap ex- 
ponent. For example, this would imply that the particle 
density of ID and 2D bosonic gases behaves as 



p(x;t) » l^ M V N {xl^,tl^,l f /l ) 



(69) 



This scaling Ansatz will be confirmed by the time evo- 
lution of a ID gas of impenetrable bosons, see Sec. lVIIB] 
We also expect that, like the equilibrium behavior, ID 
impenetrables boson gases and ID BH models of N par- 
ticles share the same TSS of the off-equilibrium dynamics 
after instantaneous variations of the trap. 



B. Power-law time dependence of the confining 
potential 

A non trivial time dependence of the confining poten- 
tial makes the issue more complicated. In the follow- 
ing, we consider a power-law time dependence, such as 
that given by Eqs. (0 and ((6]). We are again interested 
in the asymptotic trap-size dependence of the quantum 
time evolution for large initial trap size Iq. 

We use RG scaling arguments to infer the scaling be- 
havior of the time dependence of generic observables un- 
der a change of the confining potential. For this pur- 
pose, we write the perturbation associated with the time- 
dependent confining potential to the continuum theory 



i.e., 



drd d xvP T q \x\ v \<\>{x,T)^ 



(70) 



where r indicates a time variable. We are interested in 
the scaling behavior at fixed N, large Iq, with N/lo <C 1. 
The RG arguments of Sec. HJ see also Ref. [25j], may be 
extended to allow for the presence of a time-dependent 
perturbation (I70|) . and derive an off-equilibrium scaling 
Ansatz. A standard analysis of the RG dimensions of 
the coupling u leads to [56[ y u — (2 + p + zq)/p. It 
is convenient to introduce the initial trap size at t = 0, 
Iq = 1/u, with the corresponding RG dimension 



P 



2+p + zq 



(71) 



Let us consider an operator O whose low-density criti- 
cal behavior of its matrix elements is described by the RG 
dimension y in the homogeneous system. In the pres- 
ence of a chemical potential, the simplest Ansatz for the 
large-^o off-equilibrium behavior, which may be derived 
from the above RG scaling arguments, is 



The corresponding Ansatz for the low-density TSS at 
fixed particle number N is 

(0) N (x;t) a l-^ eo O N (xlo e °,Tl^ e °). (73) 

For example, the application to the one-particle density 
matrix reads 

Pi( Xl ,x 2 ;t) « lo e °M N (x i lo eo ,Tlo* eo ). (74) 

We warn that these scaling behaviors neglect pos- 
sible relevant effects related to the initial conditions, 
which may not allow us to take the Iq — > oo limit of 
the scaling functions On, more precisely of the product 
If," 6 " (O) i>f(x;t) after the variable rescalings X = xIq B ° 
and Z ee Tl~ ze °. 

The above scaling Ansatz can be reexpressed in terms 
of the instantaneous trap size 

l(t) = lor-"/?. (75) 

Replacing it in Eq. (|73p . we can write 

(0) N (x;t) ~ l(t)-y° 6 d N (xl(t)- ,Tl(t)- ze ). (76) 



The RG arguments leading to the scaling Ansatz ([72]) 
and (|73[) are qui te g eneral and can be applied to other 
models. In Ref. [56j they were applied to the XY chain 
in a space- and time-dependent trasverse field. In the 
following we challenge them against the off-equilibrium 
evolution of ID bosonic particle systems. 



VI. N PARTICLES IN A SLOWLY 
TIME-DEPENDENT TRAP 

We here discuss the behavior of N particles in a time- 
dependent confining potential, which varies slowly, i.e., 
with a large parameter t q in Eq. (|S|), and for sufficiently 
large trap sizes to be in the low-density regime. More 
precisely, we assume that the external potential is slowly 
varied so that the trap size slowly increases, correspond- 
ing to the limit t q — > — oo in Eq. ©, thus l(t) — > oo for 
t^\t q \. 



A. Adiabatic evolution 

In the case of slow changes of the Hamiltonian parame- 
ters, the system undergoes a quasi-equilibrium dynamics, 
i.e., starting from the ground state at t = 0, the evolution 
of the system passes through the instantaneous ground 
states of the BH Hamiltonian with the confining poten- 
tial V(r,t) and trap size l(t). We write the solution of 
the Schrodinger equation, 



(0)(n,x;t) « Z Q -^M (VV^/4* °)- (72) 



(77) 
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in terms of the instantaneous eigenstates (j> n of the time- 
dependent Hamiltonian (whose spectrum is discrete for 
any finite trap size), where 4> n (t) are solutions of 



H(t)<p n {t) = E n (t)cf> n (t). 



(78) 



Starting at t = from the ground state of the Hamilto- 
nian at t = 0, i.e., ^(0) = <fio(0), and writing 



l(i) = c- ie »(^a»(i)^(t) 



(79) 



where 



9„(t)= / E„{t')dt', (80) 
Jo 

the zero-order adiabatic approximation gives 

a n (t) = S n0 . (81) 

Note that the adiabatic quasi-equilibrium evolution re- 
quires the absence of degeneracies and level crossings 
during the process, but the instantaneous gap vanishes 
when t = 1 + t/t q — > 0. Thus we may already expect 
that, approaching the time corresponding to t = 0, the 
adiabiatic condition breaks down at some point of the 
evolution. We return to this point later. 

Under the quasi-equilibrium dynamics due to the slow 
increasing of the trap size, the particle number N is con- 
served because the particle number operator commutes 
with the time-dependent Hamiltonian. Since the system 
passes through equilibrium ground states, we can use the 
results obtained for the equilibrium TSS, see Sees. HI1 Hill 
and llVl The adiabatic evolution of a generic observable O 
can be obtained by computing its expectation values over 
the instantaneous ground states. After a sufficiently large 
time, when N/l(t) <C 1, we are in the low-density regime, 
thus the adiabatic time-dependence of the observables is 
obtained from the static low-density behaviors, such as 
(|43|) and ([45]) , by replacing I with the instantaneous trap 
size (|75|) . i.e., 



(O) adiab (x;0 ~ i(t)-y° 9 o N (xi(t)- e ), 



(82) 



whee l(t) is the instantaneous trap size ([75]) . Note that 
this is compatible with the dynamic TSS derived in 
Sec. El cf. Eqs. ([76]) and ([73]). 

At large N, we may use the relation ([27]) to define a 
time-dependent chemical potential fj,(t) at any t, by re- 
placing / with l(t), along the quasi-equilibrium evolution. 
For example in the ID HC-BH model, since l(t) — > oo for 
t — > oo and therefore N/l(i) — > 0, we have that fi(t) — > 1, 
which is the location of the low-density to empty-state 
transition. Within the adiabatic dynamics, the time be- 
havior of the observables related to the ground state can 
be read from that at equilibrium, by replacing the in- 
stantaneous trap size l(t) and chemical potential /x(t), 
obtained from Eq. (|27]) . in the corresponding TSS for- 
mulae obtained in Ref. (55[. It is then convenient to 
define 



which tends to zero from below in the large-t limit (after 
t q — » — oo). Asymptotically, when \p,(t)\ -C 1, the time 
dependence corresponds to varying the trap size l(t) so 
that 

\mW? e = bN 29 {l + 0[(JV/Z) 28 >]} , (84) 

where b is a p-dependent constant which can be easily 
derived from Eq. (|27p. for example b — 1 for p = 2. Note 
that the l.h.s. of Eq. (f84|) corresponds to the rescaled 
chemical potential /j, r = l 26 fi, cf. Eq. ([3"TT) . and that it 
remains constant during the adiabiatic changes since the 
r.h.s. remains fixed, apart from suppressed corrections. 



B. First-order adiabatic perturbation theory and 
breaking of the adiabatic condition 

We may also consider the first-order correction to 
Eq. (j8ip within the adiabatic perturbation theory, see, 
e.g., Refs. [H, [SSi- The first-order approximation of the 
coefficients for n > of the expansion (fT^j) over instan- 
taneous bases is 

a n {t) « -e iAe »°« f dt 1 (n\d t '\0)e~ l ^ o{t '\ (85) 
Jo 

where AO nm = 0„ — m , and, assuming nondegenerate 
states, 



(n\d t \0) 



(n\d t n\0) 

E n {t)-E (ty 



(86) 



In the low-density regime, we can use the equilib- 
rium TSS developed in the previous sections, to evaluate 
the first-order adiabatic approximation of the coefficients 
ct n (t). Energy differences behave as 



E n (t) - E (t) = ej(ty 



(87) 



where e„ generally depends on the instantaneous eigen- 
state |n). The scaling behavior of the matrix element 
(n\dtH\0) is computed considering t as a parameter. We 
evaluate the matrix element between the ground state |0) 
and one of the excited states \n) in the low-density and 
TSS limit. We generally expect 



l \d t J2[xi/l(t)] p blb i \0)=g n d t [l(ty 



(88) 



m = (i(t) - 1, 



(83) 



where g n is a (eigenstate-dependent) constant. Some of 
the lowest excited states, and in particular the lowest 
one, are obtained by exciting only the particle with the 
highest energy in the ground state, from the one-particle 
state fc = iV-lto£; = iV-l + 2j with j > 0. In the 
case of ID HC model we have 

i i 

^r ze [ dX(X p /p)<p N _ l+2j {X)y N _ x {X), (89) 
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where N is the number of particles, and the functions 
<Pfc(X) are the solutions of Eq. ([M]). Thus the integral is 
finite. It increases as N 26 at large N. In particular, for 
p = 2 

dX{X 2 /2)y N+ i{X)y N _ x {X) = j[l + OiN- 1 )}. 

The matrix element (e\dtH\0) is then obtained by per- 
foming the time derivative of the r.h.s. of Eq. in 
agreement with Eq. These results apply also to a 

ID gas of impenetrable bosons. 

Inserting Eq. (157)) and Eq. (|5H|) in the first-order adia- 
batic expansion (|85[) of the coefficients a„(i), and defin- 
ing the scaling variable Z = tIq zS °, where 9q is the off- 
equilibrium trap exponent (|7T|) . we obtain 



,t q Z D+L /b / 



d( bg n 



-it q e„C + /b 



o C e„ 



(91) 



where b = zqO /p. This expression agrees with the scaling 
Ansatz (f72|l . 

Note that it diverges logarithmically when Z — > 0, 

a n (t) « ^ ln(Z/Z ) ~ ln[/(t)/l ] ~ lnr. (92) 

Since the adiabatic perturbative expansion requires 
|cKn(t)| <C 1, it fails when r becomes too small. This 
is not unexpected because when r — ¥ the spectrum 
tends to be degenerate. 

A simple example of the breaking of the adiabatic evo- 
lution when approaching a Hamiltonian with vanishing 
instantaneous gap is provided by a quantum oscillator 
with a time-dependent frequency, see App. [Cj 



VII. ID IMPENETRABLE BOSONS IN A 
TIME-DEPENDENT HARMONIC TRAP 

In this section we determine the trap-size dependence 
of the off-equilibrium evolution of a ID gas of impen- 
etrable bosonic particles in a time-dependent confining 
harmonic potential, i.e., p = 2, starting from an equilib- 
rium ground state configuration with initial trap size Iq. 
We consider instantaneous changes to a confining poten- 
tial with different trap size If, and also power-law time 
dependences such as 



V(x,t) = ^K{t)x 2 , 



(93) 



where 

n{t) = KQT q 



l{tf 



T=l+t/tq, KQ = l/ll (94) 



and t q is a time rate. In the following analyses of the 
power-law time dependence, we set t q = 1 for simplicity. 



A. Off-equilibrium time evolution 

As shown in Ref . [3l| , see also [4!| , the time-dependent 
wave function of the system can be derived from the so- 
lutions tjjj (x, t) of the one-particle Schrodinger equation 



id t ipj(x,t) 



1 



d 2 x + V(x,t) 



(95) 



with the initial condition tpj(x, 0) = 4>j(x) where <fij(x) 
are the eigensolutions of the Hamiltonian at t = 0, 
characterized by a trap size Iq, with eigenvalue Ej = 
(j + 1/2) ci. Eqs. ([6Tj). The solution can be ob- 
tained introducing a time-dependent function s(t), writ- 



ijjj(x,t) = s 1/2 (t>j(x/s) x 



(96) 



. sxr 

xexp | % ihiA 

5 2s 3 



l dt' 



where 4>j( x ) is the j eigenfunction of the Schrodinger 
equation of the one-particle Hamiltonian at t = 0, thus 
with trap size Iq, and s{t) satisfies the nonlinear differen- 
tial equation 



s + n(t)s — kqs 



-3 



(97) 



with initial conditions s(0) = 1 and s(0) = 0. 

The time-dependent wave function $ of N impenetra- 
ble bosons, with §(x, 0) given by the ground state of the 
Hamiltonian at t — 0, can be obtained following the same 
steps as at equilibrium, see Sec. IIV) obtaining 32] 

$(xi, ...,x N ;t) = A{xi, ...,x N )^{xi, ...,x N ;t), 
V(xi,...,x N ;t) = -^det^ifo;*)], (98) 



where the determinant involves the N lowest eigensolu- 
tion a fixed t. Then, using Eq. (|96[) . one can write the 
wave function of the ground state of an iV-particle system 
as 

<^{x 1 ,...,x N ;t) = s- N/2 $(x 1 /s,...,x N /s;0) x 



(99) 



where <&(x\, xjv; 0) is the wave function of the ground 
state for the Hamiltonian at t — 0. The time-dependent 
one-particle density matrix reads [4{| 

Pi(x,y;t)= (100) 
= N j $(x,x 2 ,.:,x N ;t)*$(y,X2,...,x N ;t)dx 2 ...dx N 



*f(</ 2 -- 2 ) 

Is 



' I pi(x/s,y/s;0)exp 



where pi(x,y;Q) is the equilibrium one-particle density 
matrix at a trap size I = Iq, i.e., 



pi(x,y;0) = pi(x,y)\ l=lo 



(101) 
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s(t) 




B. TSS at instantaneous quenches 

We now show that the time evolution after instan- 
taneous changes of the trap size is consistent with the 
scaling Ansatz put forward in Sec. IV Al in terms of the 
equilibrium trap exponent 8 = 1/2. 

Let us first consider an instantaneous drop of the trap. 
The energy after the quench can be computed within ID 
BH model in the TSS limit, by evaluating the expectation 
value of the unconfined BH Hamiltonian H u over the 
ground state |0 C ) of the confined BH Hamiltonian 



H c =U u + S ^V{x l )n l 



(107) 



FIG. 14: (Color online) The function s(t) for kq = 1, cf. 
Eq. GH). 



computed in Sec. El cf. Eq. (|55|) . 

Examples of explicit solutions of the function s(t), cf. 
Eq. (j9"7| . are the following. 

(i) Instantaneous drop of the trap, so that n(t) = for 
t > 0, 



?(t) = v 7 ! + Kot 2 



(102) 



(ii) Instantaneous change to a confining potential with 



trap size If, so that n(t) — l f for t > 0, 



s(t) = J 1 + (r 2 - 1) sin(Kg /2 t/r) 



(103) 



where r = If/lo- 

(iii) Linear time dependence of the trapping potential, 
i.e., q = 1 in Eq. ([M|. 



s(t) = [ReW(r)] 



where the complex function W(t) is the solution of the 
differential equation 



iW' = nl /2 {W 2 -t) 



(105) 



with W(l) — 1, 5 which can be written as a combination 
of Airy functions, [52] 



W(t) = in- 



„i/ 6 Bi'(-KQ /3 T) + cAi(-kI /3 t) 



Bi(-Ko /3 r) + cM(-k /3 t) 
^ /6 Ai(-^ /3 )- ? Ai'(-^ /3 )- 



(106) 



A plot of s(t) is shown in Fig. H4l 



5 In the general case, i.e., when V(x;t) = k(£)x 2 /2, the replace- 
ment Ill04t leads to the differential equation iW 1 = reJ^W 2 — 
k~ 1/2 k(t) with W(l) = 1. 



in the low-density region N/l -C 1. We have 



^ = (0 C |H„|0 C ) 
For p = 2, we have 



7(si)ni|0 c ). (108) 



JV-l 



Ei — I, 



o 1 [E( fc + 1 /2) 



fc=0 



N-l 



dA(A 2 /2)> w k (X) 2 =-n-. (109) 



fc=0 



TV 2 
4^ 



For generic values of p, we have Ei ~ jyio+i jpe xhcrc- 
fore, for large initial trap size Iq, thus -/V/io <C 1, only 
low-energy states are involved. 

We again expect that in the low-density regime the 
asymptotic trap-size dependence is that of the gas of im- 
penetrable bosons, and that the lattice structure of the 
BH model gives only rise to suppressed power-law cor- 
rections. Therefore, in the case of the harmonic trap we 
can use the general solutions reported in the previous 
subsection to derive the TSS behavior at a quench. 

The time-dependence of the one-particle density ma- 
trix, after turning the trap off, is obtained by inserting 
the function s(t) of Eq. (pXJ2|) into Eq. ([TOO)) . Then, using 
equilibrium relation 



Pl (x,y;0) n l- 9 M(x/l e ,y/l 9 ), 



and defining 



X = x/ll Y = y/ll Z = t/F 6 



(110) 

(111) 
(112) 



Q(Z) = Vl + Z 2 , 
where 9 = 1/2 is the equilibrium trap exponent, we write 



pi(x,y;t) = l e Q 1 x 
xM N {X/Q,Y/Q)exp 

The particle density is given by 



(113) 



X 2 ) 



p{x- 1) = Pl {x, x; t) = l^Q^VNiX/Q), (114) 
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FIG. 15: (Color online) Iqp(x; t) for some values of Z = Iq t, 
for N — 10 (below) and in the limit N — s- oo (above), in the 
case of a quench to the unconfined Hamiltonian. 

where T>n(X) can be derived from Eqs. (|40|) and (|46|. 
Analogously, one can derive the particle-density correla- 
tion G n , cf. Eq. (jtJo) . obtaining 

G n (x,y;t) = lo 2f> Q- 2 g N {X/Q,Y/Q). (115) 

Results for the particle density and the one-particle den- 
sity matrix are shown in Figs. [Pol and [Pol 

In the case of a quench to a larger trap size If > Iq, 
replacing Eq. (|103|) into Eq. (|100p , we again obtain the 
expression (|113p . but 

Q{Z) = ^l + (r 2 -l)[sin(Z/r)] 2 , (116) 

where r = If/lo- Therefore, we have a periodic time 
evolution with period Z p = nr. Fig. [T7] shows results 
for the periodic time evolution of the particle density for 
N = 10 and r = 2. 

For a large number of particles in a harmonic potential, 
we can derive the TSS using the asymptotic behavior 
given by Eqs. flBl]) and (|B2]) . We obtain 

^-h^f-WQr (117) 




FIG. 16: (Color online) Iq \pi (0, x; t)\ for some values of Z = 
l^ 20 t, for iV = 10 and in the case of a quench to the unconfined 
Hamiltonian. 




FIG. 17: (Color online) Time dependence of the rescaled 
particle density lap(x;t) for some values of Z = lQ 29 t, for 
N — 10 and in the case of a quench to a trap with size I / = 2/o • 
It oscillates between the Z = and Z — ir curves. 

The above results show that, after instantaneous 
changes of the trap size of the harmonic confining poten- 
tial, the trap-size dependence satisfies the scaling Ansatz 
put forward in Sec. IV Al 

It is worth noting that the time-dependence of the par- 
ticle density p(x;t) and its correlation G n (x,y;t) can be 
reexpressed as their equilibrium TSS with an effective 
time-dependent trap size 

l(t) = l aS {t) l ' e = loQ(Z) 1 / 6 , (118) 

so that 

p(x;t) = l{t)- e V N {X), X = x/l(t) e . (119) 

We finally mention that the case of a gas of impene- 
trable bosons in a hard-wall trap, and its expansion after 
the drop of the trap, was considered in Refs [7l|, [72[. A 
hard-wall trap of size L corresponds to the p — > oo limit 
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of the confining potential, cf. Eq. ([2]), with trap size 
I = L/2. One can easily check that the time evolution of 
the particle density after the drop of the trap, computed 
in Ref. [72|], is consistent with the scaling Ansatz 
taking into account that the p — > oo limit of the trap 
exponent (I10|l is 8 = 1 and Iq = L/2. 



C. Power-law time dependence of the trapping 
potential 

We now consider the case of a power-law time depen- 
dence of the confining potential, cf. Eq. 
Let us define the quantities 



S(Z) = lf o/2 s(t), 



= 7- 2e o 

<T1 



(120) 



Q. 




where 



FIG. 18: (Color online) Time dependence of the particle 
density for N — 10 in a time-dependent trap with «(t) = r, 
corresponding to lo = 1. 



<l 



(121) 



is the off-equilibrium trap exponent obtained by replacing 
z = 2 and p = 2 in Eq. (|TTj) . S(Z) satisfies the equation 



S" + z q s = s- 



(122) 



where 5"' = d\S. Then, using the equilibrium relation 



([TTfJ]) . we rewrite Eq. ([TOO)) as 

p 1 (x,y;t) = lo 8o S- 1 x 
xMn{X/S,Y/S) x exp 



(123) 



i—(Y 2 
2S v 



where X = x/l ° , Y = y/l ° 7 and .Mat is the same scaling 
function appearing in Eq. (|67j) . 

The evolution of the particle density is easily obtained: 



p(x; t) = Pl (x, x- 1) = lQ e °S- l V N (X/S), 



(124) 



where Pjv is the scaling function (|46l) . For a large number 
of particles we can derive the off-equilibrium TSS using 
the asymptotic behavior given by Eqs. (jBlj) and (|B2[) . 
We obtain 



(27V) 



1/2 



1 



X 2 



2AS 2 ' 



(125) 



The time dependence of p(x\ t) can be again reex- 
pressed using the equilibrium expression with an effective 
trap size 



f(t) = f^s 1 ' 6 = i a -s(t) 1/e . 



(126) 



The function l(t) increases monotonically with decreasing 
t < 0. If one prefers to invert the time evolution, so that 
the effective trap size increases with increasing t > 0, it is 
sufficient to redefine r = 1-t in Eqs. (JM]) and (TIM)) . Note 
that l(t) remains finite for r = 0, i.e., when the external 
potential (|93l) vanishes, indeed s(r = 0) ~ 1.14313, then 



it diverges in the limit t — > — oo (note that for r < and 
g = 1 the potential changes sign, so it does not trap 
the particles anymore). 

The above scaling behaviors are apparently consistent 
with those predicted by the scaling arguments of Sec. lVBI 
for the off-equilibrium TSS in the low-density regime. 
However, the function S(Z) maintains a residual depen- 
dence on Iq, beside on Z, due to the initial condition of 
s(t) which corresponds to l~ q6o/2 S^ 29 ") = 1. Thus, the 
scaling Ansatz can be actually considered as fully verified 
only if the function S(Z) has a nontrivial scaling limit for 

Iq — > OO. 

In the case of a linear time-dependence of n{t), i.e., 
n(t) = kqt, the function S{Z) can be derived from the 
corresponding solution s(t), cf. Eq. (|104p . Then, using 
the equilibrium results of Sec. IIV1 we obtain the time de- 
pendence of the one-particle density matrix, the particle 
density, particle density correlators, momentum distribu- 
tion, etc.... Some results for the particle density and one- 
particle density matrix are respectively shown in Figs. 1181 
and[H 

An important remark is in order. The analytical solu- 
tion in the case of a linear time dependence shows that 
the function S(Z) does not have a nontrivial scaling limit 
for Iq — > oo, indeed it appears to diverge, roughly as 
(InZo) 2 at fixed Z. This may reflect the fact that the 
initial conditions are somehow weakly relevant, leaving 
a residual weak (logarithmic) dependence in the large- 
Iq limit. We should further note that S(Z) = s(t) with 



Z = L 



-20, 



t (with r = 1 +t) has a nontrivial Iq 



limit 



satisfying the differential equation S" + ZS = 0. But this 
rescaling would not fit any scaling behavior consistent 
with the dynamic exponent z — 2. This point deserves 
further investigation. 
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FIG. 19: (Color online) Time dependence of the absolute 
value of the one-particle density matrix for N = 10 in a time- 
dependent trap with n(t) = r, corresponding to Iq = 1. 



VIII. SUMMARY AND CONCLUSIONS 

We study the trap-size dependence of the quantum 
behavior of dilute gases of bosonic particles in the 
presence of a confining potential trapping the particles 
within a limited spatial region. We consider systems of 
bosonic particles constrained in an optical lattice, de- 
scribed by the Bose-Hubbard (BH) model in the pres- 
ence of a confining potential coupled to the particle den- 
sity, cf. Eq. ([T]). In the case of a harmonic potential 
V{x) — v 2 x 2 /2, the corresponding trap size is defined as 
I = \fj jv where J is the hopping parameter. We con- 
sider systems at equilibrium and off equilibrium during 
the unitary time evolution arising from changes of the 
trapping potential, at zero temperature, i.e., at a suffi- 
cient low temperature to neglect its effects. We investi- 
gate the trap-size dependence in the low-density regime 
using the framework of the trap-size scaling (TSS) the- 
ory til nu. 

Using scaling arguments, we infer the power-law trap- 
size dependence of observables related to the equilibrium 
lowest states of a dilute gas of N particles in the low- 
density regime. The low-density regime of the BH model 
dJ), i.e., Na d /l d « 1 where a is the lattice spacing, 
can be seen as the critical regime of a quantum tran- 
sition from low density to the empty state, which may 
be considered as a n — Mott transition. Mott transi- 
tions driven by the chemical potential are described by 
the nonrelativistic $ 4 continuum theory where the 
dynamic exponent is z = 2 and the RG dimension of 
the chemical potential is y M = 2, in one and two spa- 
tial dimensions. In the presence of a confining potential, 
the power-law trap-size dependence is described by the 
equilibrium trap exponent [H [55[ 9 = p/(p + 2) where 
p is the power of the confining potential. This allows us 
to derive the universal scaling features of the asymptotic 
power-law trap-size dependence keeping fixed the particle 
number N. For a generic observable, whose low-density 



critical behavior is described by the RG dimension y Q in 
the homogeneous system, we obtain the scaling Ansatz 
(0) N (l,x) « l-*° O N {xl-% see Sec. ED 

The equilibrium TSS scenario is verified in ID systems 
by analytical and numerical calculations. We show an- 
alytically that the expected TSS holds in the hard-core 
(HC) limit U oo of the BH model. We compute the 
scaling functions of some observables, such as the particle 
density and its correlators, the one-particle density ma- 
trix, see Sec. IIII B 21 The universality of the low-density 
TSS with respect to the on-site repulsion coupling U is 
supported by numerical calculations at a finite value of 
U, i.e., (7 = 2, using DMRG methods. We show that 
the asymptotic TSS of N particles at equilibrium de- 
scribed by the ID BH model, in the HC limit and at 
finite U, is identical to that of a ID gas of impenetrable 
bosons (Tonks-Girardeau model), with appropriate def- 
initions of the trap size (in the case of harmonic traps 
the trap sizes are proportional to the inverse frequency 
in both models). The lattice structure gives rise to sub- 
leading 0(l~ ) scaling corrections in the HC limit of the 
BH model. The approach to the asymptotic behavior is 
slower at finite values of U, in agreement with the RG 
arguments which predict subleading 0(l~ e ) scaling cor- 
rections. We argue that the same scenario applies to the 
Lieb-Liniger model with a finite contact interaction in 
the low-density regime, i.e., it presents the same univer- 
sal asymptotic TSS with 0(1 ) scaling corrections. 

We investigate the trap-size dependence of the off- 
equilibrium dynamics due to time-dependent confining 
potentials, such as V(r,t) ~ (1 + t/t q ) q r p , or instan- 
taneous changes of the trap size, including the in- 
stantaneous drop of the trap. We extend the scaling 
Ansatz for the trap-size dependence at equilibrium to 
off-equilibrium quantum evolutions, to describe the TSS 
with respect to the initial trap size Iq, see Sec. [V] We 
argue that (0) N (x;t) « l^ Vof> O n {xIq 6 ,tl^ z6 ,l s /l ) in 
the case of an instantaneous change of the trap size 
from Iq to If . In the case of a power-law time depen- 
dence, we introduce an off-equilibrium trap exponent, 
given by 8 = 1/(2 + q) in the case of a harmonic 
trapping potential, and put forward the scaling Ansatz 
(0) N {x;t) « l- v ° ea O N {xlQ eo ,Tl- z00 ). The above results 
are expected to be quite general in the dilute regime of 
ID bosonic gases, such the lattice BH model and contin- 
uous Lieb-Liniger model. 

We then analyze the trap-size dependence of the off- 
equilibrium dynamics of ID bosonic gases with respect 
to the initial trap size Iq, using adiabatic approximations 
in the case of slow changes of the parameters, and ex- 
act solutions of the Schrodinger equation of N impene- 
trable bosons in time-dependent traps or after instanta- 
neous changes of the trap size. The evolution after in- 
stantaneous quenches agrees with the corresponding off- 
equilibrium Ansatz (|68|) . where the equilibrium trap ex- 
ponent 9 characterizes the power-law dependence on the 
initial trap size, see Sec. lVII 5] In the case of a power-law 
time dependence of the potential, the evolution supports 
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the scaling Ansatz in terms of the off-equilibrium trap 
exponent 9q, see Sec. I VII CI However, in the case of a 
linear time dependence, for which we have an analyti- 
cal solution, the simplest Ansatz ([75)) does not provide a 
complete description of the asymptotic large-^o behavior 
because the resulting scaling functions maintain a weak 
logarithmic dependence on Iq in the large-/o limit, demon- 
strating that the initial conditions are somehow relevant. 
This point deserves further investigation. 

Our results are of experimental relevance for systems 
of cold atomic gases trapped by a confining potential. 
Indeed, the easy tunability and long characteristic time 
scales of these systems may allow a careful study of the 
trap-size dependence of the zero-temperature properties 
of iV-particle boson gases, in the continuum and on opti- 
cal lattices, at equilibrium and off equilibrium in a time- 
dependent confining potential. 

Helpful discussions with P. Calabrese, D. Giuliano, M. 
Mintchev and G. Morchio are gratefully acknowledged. 



Appendix A: TSS in traps induced by a spatial 
dependence of the hopping parameter 

As suggested in Ref. [73]], ultracold atomic systems in 
optical lattices may be also get trapped by appropriate 
spatially inhomogeneous hopping parameters in the BH 
model. An example is given by the model 

Uu< = ~ £ + 4bj) + % £ Mm 1) (Al) 



0.5 



2^2 



with 



h{x) = 



X i I X j 



p 



-(x/iy 



2 

-l 



(A2) 



and Xi are the positions of the sites of the lattice. The 
rescaled hopping parameter tij tends to one at the middle 
of the trap and vanishes at large distance, giving rise to 
an effective trap, with trap size I. 

In the HC U oo limit, the Hamiltonian can be diag- 
onalized exploiting the fermion quadratic representation 
U23J) with 



(A3) 



following the procedure outlined in Sec. MI B 11 cf. 
Eqs. ©2), 433) and 43D). 

In the dilute region, i.e., for sufficiently small N/l, we 
can follow the same steps of Sec. lIII B~T1 to arrive at a con- 
tinuum TSS limit. We end up with the same Schrodinger- 
like equation (|39|) after the same rescalings (|36l) , (|37|) and 
(08) and = p/{p+2) as well. Thus, the TSS arising from 
the spatial-inhomogeneity of the hopping parameter, like 
Eq. (|A2I) . is identical to that of model (TTJ), i.e., of a trap 
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FIG. 20: (Color online) The large- N behavior of V N (X) for 
p = 2. The full line shows its N -> oo limit, cf. Eq. (jM) 
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FIG. 21: (Color online) Comparisons of V N (X) for p = 2 
and iV = 10 and N = 40 with their next-to-leading approxi- 
mations, cf. Eq. (|B3[) . denoted by "1/N" in the figure. 



achieved by coupling an external potential to the particle 
density. One can also infer that scaling corrections are 
0(l- 26 ) as well. 



Appendix B: Trap-size scaling for a large number of 
particles 

In this appendix we determine the largc-A behavior 
of the TSS functions of the observables considered in 
Sec. HUB 21 



1. The large- iV behavior of the TSS functions 

In the case of the harmonic potential, the rescaled par- 
ticle density, cf. Eq. (14^1) . behaves as 



V N {X) = N 1 ' 2 



r d{ x) + ^1 + o{i/n 2 ) 



,(B1) 
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FIG. 22: (Color online) The scaling functions T> N (X), cf. 
Eq. (I45p . for p = 4 (below) and p — » oo (above). 



where X = N 1 ' 2 X. The approach to the large-iV be- 
havior is shown in Fig.[2Hl where N~ 1 / 2 Djq{X) is plotted 
versus X for several values of N . The leading \axge-N 
behavior is 



R D (x) = -y/2 



(B2) 



for x < \/2, and Rd(x) = for x > \2. This is also 
the N — > oo limit of the particle density in a bosonic gas 
of impenetrable bosons [4l], H3] ■ Actually, since the low- 
density limit of the HC-BH model matches the behavior 
of a gas of impenetrable bosons, as discussed in Sec. IIV| 
we can use results obtained for the TG gas [4l| to infer 
that 



(- l^cos [iVg(a;)] 



C D {x) = 

q(x) — x\J 2 — x 2 + 2arcsin(x/v / 2). 



tt\/2(2 



(B3) 



A comparison of T>n(X) for N = 10 and N = 40 with 
their next-to-leading large- N approximations is shown in 

Fig. 03] 

Analogous results can be derived for other power laws 
of the confining potential. Results for p = 4 and p — > oo 
are shown in Fig. [5H The particle-density scaling func- 
tions show again N peaks, with an underlying structure 




FIG. 23: (Color online) Plot of N^GnIX, Y) for X = Y/2, 
X = 0, and X = —Y, vs N 1/2 Y, for several values of N. 
They approach unique curves with increasing N. Even and 
odd values of converge from opposite sides. 



scaling as V N {X) w iV 2 / a i? D (X/iV 1 / 3 ) for p = 4 and 
T>n(X) = N/2 + 0(1) for p —> oo (at least not too close 
to X = 1), with the oscillatory terms suppressed by 1/N 
with respect to the leading terms. These results suggest 
the general behavior 



V N (X) = N e [RoiX/N 1 - 6 ) + 0(l/N)] 



(B4) 



for any power p, where the function Ru depends on p. 

Fig. [53] shows results for the scaling function Qn(X, Y) 
associated with the density-particle correlation, cf. 
Eq. (fTTj) . for several values of N, obtained from Eq. (|48|) 
for the harmonic potential. These plots show that the 
large- iV behavior is 



G N {X,Y) « NR G {N 1/2 X,N 1/2 Y) 



(B5) 



Note the different 7V-rescaling of the spatial coordinates 
with respect to that of the particle density. 

Concerning the one-particle density matrix, cf. 
Eq. (fl~3|) . we note that the scaling behaviors (|B4j) and 
(|52p. of the particle density and the momentum distri- 
bution respectively, can be both derived from the follow- 
ing nontrivial large- N scaling behavior of the one-particle 
density matrix: 

Pl (x, y) « (N/l) 9 B[N- 1+e X, N\Y - X)}, (B6) 

where B is a scaling function, X = x/l and Y = y/l 
(note the different power of N in the two arguments of 
the function B; in the case of a harmonic potential they 
are N^^X and 7V 1 / 2 (F - X) respectively). The above 
scaling behavior would imply that, with increasing N, 
the region where the diagonal component is significantly 
nonzero increases as N 1_e , while the width around it de- 
creases as N~ 6 . The approach to this large- N limit is 
generally characterized by 0{N~ 1 / 2 ) oscillating correc- 
tions in the case of the harmonic potential. 
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2. Power-law behavior at the boundaries where 
the particle-density vanishes 

In this section we discuss the asymptotic large-TV 
power-law behavior at the boundaries of the trap where 
the particle density is suppressed. 

The large-TV asymptotic behavior found in App. 151 
holds for X = N~ 1/2 X < X c where X c = ^2 is the value 
where Rjj(X) = 0, which is the point around which the 
particle density vanishes in the large- N limit. However, 
around the spatial points corresponding to X c another 
power law behavior sets, as suggested by the behavior 
around X c of the curves shown in Fig. [21] This fact was 
already noted within the Gaussian unitary ensembles of 
random matrices [50l r74j , whose eigenvalue density cor- 
responds to the particle density in harmonically trapped 
systems of impenetrable bosons. An analogous change of 
power law is observed at fixed chemical potential, thus 
N ~ I, at the boundaries of the trap [55| • 

This phenomenon is related to a real-space transition 
between the low-density particle regime, for X < X c , 
and the empty state for X > X c , which occurs at the 
points x c corresponding to X c . Thus we expect that the 
region around x = x c develops critical modes related to 
a low-density Mott transition. The effective external po- 
tential at x c can be obtained by expanding the trapping 
potential around x c , thus obtaining an approximately lin- 
ear potential Vi(x) ~ x — x c . Around x c , other critical 
modes develop with length scale £ ~ where a is the 
exponent associated with a linear external potential. The 
value of a can be inferred by RG arguments analogous to 
those leading to the determination of the trap exponent 
6 at the low-density Mott transition [25|, [55| , which give 
(j = 1/3. 6 For example, this implies that 

X c - X max ~ N- 2 / 3 (B7) 

where A max corresponds to the abscissa of the rightmost 
maximum of T>n(X). More generally, we have the scaling 
behavior 

HuiN^N-^VNiN^iXc + N- 2 ' 3 z)\ = f(z). (B8) 

The scaling function f(z) can be obtained from related 
computations within the Gaussian unitary ensembles of 
random matrices [5(| [zl] : 

f(z) = 2 1 / 2 |Ai'(2 1 / 2 z)| 2 - 2z|Ai(2 1 / 2 z)| 2 . (B9) 

Fig. [M] shows that the above asymptotic behavior is 
rapidly approached in the large- N limit. 



The exponent a can be determined by a RG analysis of the per- 
turbation corresponding to a linear potential Vi(x) = ux, i.e., 
Jd d xdt Vi(x)\ij>(x)\ 2 , at the fixed point of the continuous theory 
describing the Mott transition [2111 . The exponent a is related to 
the RG dimension y u of the parameter u, which can be obtained 
from the relations y u — 1 = d + z — yuia = Vft = 2, thus y u = 3, 
and therefore a = l/y u = 1/3 for d = 1 and d = 2. 




FIG. 24: (Color online) Plot of N~ 1/6 V N [N 1/2 (X c + 
N~ 2/3 z)] versus z for p = 2, where X c = 2 1/2 , for TV = 
10, 20, 100. The curves appear to converge toward the 
TV -> oo limit ([B9]> . 



Appendix C: The quantum oscillator with a 
time-dependent frequency 

Let us consider a quantum oscillator described by the 
Hamiltonian 

^ = y + ^)y (CI) 
with a time dependent frequency 

K{t) =L0 2 ) {l+t/t q ) q =ulr\ (C2) 

where t q is the time rate of the time dependence. In the 
following we set t q — 1 for simplicity; its dependence can 
be easily inferred by appropriate rescalings of the results. 
We assume that at t = the oscillator is in its ground 
state, i.e., its wave function is 

M*) = K/^) 1/4 e^ 0l2/2 . (C3) 

The evolution equation 

id t ip(x,t) = Hip(x,t) (C4) 

with ip(x, 0) = ifjoix) preserves the Gaussian spatial de- 
pendence. We write the solution of Eq. (IC4[) as 

Tp{x,t) = ( Wo A) 1/4 e- w(t)woa;2/2+ " (t) , (C5) 

where w(t) and z(t) satisfy the equations 

iw = uj (w 2 -T q ), (C6) 
iz = uiqw/2, 

with initial conditions w(0) = 1 and z(0) = 0. In the 
case of a linear dependence of re, i.e., 

K(t) = l + t = T, (C7) 
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FIG. 25: (Color online) Some results for the coefficients of the 
expansion (|C9[) versus t, using exact results from Eq. (|Clip 
and from the adiabatic approximation (|C12[) . 



the solution is w(t) — W{t) where W(t) is the complex 
function given in Eq. (|106|) . and 



z(t) = 



ILOQ 



dt' w{t') 



(C8) 



Note that tp(x,t) remains exponentially suppressed at 
large x even at t = — 1 when k = 0, indeed uj(— l)|o;o=i — 
W(0)L=i = 0.765265 + zO.346358. 

It is interesting to compare the exact solution (|C5j) 
with the evolution predicted by the adiabatic perturba- 
tion theory. In the following calculations we set ljq = 1. 
We expand the wave function as 



where £ n K(*)| 2 
eigenstates, i.e., 



ij){x,t) =J2a n (t)0 n {x,t), (C9) 

n 

1 and </> n (x,t) are instantaneous 

(CIO) 



H<j) n (x,t) = E n (t)4> n (x,t), 
E n = tj{t){n + 1/2), uj(t) = 



4> n {x,t) 



1 



7r 1 /4(2"n!) 1 / 



Yj^H n [xuj(t) 



l/2i-uj(t)x' 1 /2 



where H n (x) are the Hermite polynomials. The coeffi- 
cients a n (t) are given by 



a n {t) = / dxip(x,t)*(f> n (x,t). 



(Cll) 



Note that a n — for odd values of n (when the corre- 
sponding eigenfunction (p n is odd). Moreover, a n (t) —> 
for t -> -1. 

Assuming a very slow variation of the Hamiltonian 
parameters, we expect an adiabatic evolution, i.e., the 
system starting from the ground state at t = evolves 
through the instantaneous ground states 4>o(x,t). Thus, 
the leading behavior is given by a n (t) — e ie °^^„ where 
Q n (t) = Jo E n (t)dt. The time dependent coefficients 
a n {t) can be computed to the next-to- leading order of the 
adiabatic perturbative expansion, see, e.g., Refs. |68ll69j. 
obtaining 



a n (t) 



dt' _ 
o E n (f) 



n\d t H\0) 



-iAe„ (t'] 



E (t') 



C12) 



where A© nm = 9„ — 9 m . Beside ao(t), only the coeffi- 
cient a2 (t) differs from zero in the first order adiabatic ap- 
proximation. Note that the limit t — > — 1, thus t — > 0, is 
singular; indeed | a.2 (i) | ~ | lnr|. This shows that the adi- 
abatic approximation fails approaching the critical point 
where the spectrum tends to become degenerate. 

Some results for the coefficients a n {t) are shown in 
Fig. as obtained from Eq. (|C11[) and their adia- 
batic approximation. We find that |ao(i)| 2 > 0.99 for 
t = 1 + t > 0.5, which is the region where the zero- 
order adiabatic approximation works within 1%, and 
|ao(i)| 2 + \a 2 (t)\ 2 > 0.99 for r > O.2., which is where 
the first-order adiabatic approximation is effective. 
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